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We study the dynamics of inorganic antiferroelectric nanoclnsters formed during an order-disorder 
transition and demonstrate the coexistence of the two phases in a region of 2 — 3 K around the 
transition temperature Tjv ~ 215 K. Single crystals of NH 4 H 2 ASO 4 , a model hydrogen-bonded 
compound, show an antiferroelectric-paraelectric transition studied by means of highly sensitive 
magic angle spinning NMR at 21.1 T. The phase co-existence is demonstrated by a double-peak 
structure of the chemical shift. Two-dimensional chemical exchange spectroscopy and spin-lattice 
relaxation time (Ti) measurements show that the clusters are dynamic with sizes ~50 nm and 
lifetimes approaching seconds as T —>■ Tjv. Their occupancy increases rapidly to fill the crystal 
volume below Tjv. This study provides evidence for the commonality of the phase transitions in 
systems with electric properties and provides an improved spectroscopic method for such studies. 

PACS numbers: 77.84.-s, 63.20.Ry, 76.60.-k, 77.80.B- 


An atomic-scale view of order-disorder transitions can 
significantly improve the understanding and theoretical 
modeling of interactions in the time domain, at the 
nanoscale level. Generally speaking, hydrogen-bonded 
compounds, such as the KH 2 PO 4 (KDP)-family^, are 
considered to be significant for ferroelectric (FE) tran¬ 
sitions from several points of view: (a) They serve as 
models of the phenomenon of hydrogen-bonding that is 
of fundamental importance in the physics and chemistry 
of ice and water^ as well as the structure and func¬ 
tion of biological molecules such as DNA^; (b) the KDP 
family finds extensive applications in electro-optic de¬ 
vices because of their ferroelectric behavior^; and (c) 
despite recent advances in the theoretical modeling of 
FE and related order-disorder phase transitions via com¬ 
puter simulations'^’® and analytical solutions®^®, a general 
model of the transition dynamics is not available^^’^^. We 
present a study of the antiferroelectric (AEE) phase tran¬ 
sition in an inorganic compound NH 4 H 2 ASO 4 (hence¬ 
forth ADA)^^“^® tailored to target the nanoscale and dy¬ 
namical aspects of such a transition, which is the missing 
piece of many theoretical approaches. One long-standing 
unanswered query is that while evidence of pretransi- 
tional cluster formation has been seen in the case of 
the FE transitions (in KDP systems^^^^® KH 2 ASO 4 and 
RbH 2 As 04 using ^®As NMR and nuclear quadrupole res¬ 
onance), similar measurements have failed to detect pre- 
transitional clusters during the AFE transitions^'^. Such 
a failure casts some doubt on the universality of the dy¬ 
namics accompanying the ferroelectric phenomenon in 
general, since computer models predict a similar picture 
of the phase transition dynamics in any lattice that un¬ 
dergoes an order-disorder phase transition^^®. 


The results of the above-mentioned NMR 
experiments^^^^®, together with the theoretical 
predictionssuggest that in the AD A-type AEE 
lattices, the pretransition clusters must be at the 
nanoscale, thus requiring increased spectral resolution 
and sensitivity. The ratio ry = VQjbv, where vq is the sig¬ 
nal frequency and 8v the linewidth is an accepted figure 
of merit for spectral resolution. For ADA, the previous^® 
77 was lT/(lmT)~ 10®. The recent development of 
magic angle spinning®^ (MAS) NMR techniques at 21.1 
T at the National High Magnetic Field Laboratory 
(NHMFL), gives r] -(91x10® Hz)/5 Hz-2xl0®' for 
the ®®N nucleus, constituting a four order-of-magnitude 
enhancement of the spectral resolution®®. This technique 
enables the detection of the low-symmetry AFE clusters 
(co)existing in the high-symmetry [paraelectric (PE)] 
phase of ADA. These pretransitional AEE clusters are 
dynamic and we present experimental evidence that 
atoms are moving between the clustered region and 
the rest of the lattice. The results provide support 
for theoretical models studying such dynamic atomic 
exchange and cluster formation when approaching an 
order-disorder phase transition^^®. 

We employed a ®®N-enriched ADA crystal with a size 
of 1x1x1 mm® and — 9% enrichment, prepared as re¬ 
ported earlier®®. ADA is a colorless transparent crystal 
that at room temperature crystallizes in a body-centered- 
tetragonal unit cell®®, space group lA2d, with the unique 
(c-) axis cell length = 7.709 A, and a = b =7.674 A. The 
structure changes to an orthorhombic (P 2 i 2 i 2 i) phase at 
temperatures below the AEE phase transition, measured 
to be —215 K. It becomes antiferroelectric, with the H’s 
localized asymmetrically along one of the O’s [1,10,11]. 
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FIG. 1. (Color online) (a) isotropic chemical shift of 
ADA, showing an abrupt change at Tjv =215.2 K as a func¬ 
tion of temperature. Data points indicate the maxima of 
CPMAS spectra, such as shown in (b). (b) Contour plot of 
CPMAS spectra (arbitrary units) in the vicinity of Tat, with 
darker areas indicating the location of maxima or high signals. 
The insets show three spectra; the vertical axes are of same 
range for all three and the horizontal axis is identical to the 
one of the contour plot. Tjv was defined as the temperature 
at which the PE (dark gray) and AFE (light red) peaks have 
equal integrated intensities. 


As a spectral reference, we use powdered Valine. The 
spin-lattice relaxation time Ti is measured by enhancing 
the saturation of the nuclei via a standard cross polar¬ 
ization (CP) protocob^. In addition, the MAS technique 
allows us to obtain isotropic resonance signals by spin¬ 
ning the sample around an axis tilted by an angle of 54.7° 
away from the applied magnetic field^^. At this magic an¬ 
gle, the inhomogeneous dipolar broadening is averaged 
out by the fast spinning. The NMR measurements were 
made on the ultrawide bore 900 MHz (21.1 T) NMR spec¬ 
trometer at NHMFL with a 5 kHz spinning speed and 
temperature stability of 0.1 K. In this setup, the 
and ^H Larmor frequencies are a;Ar/27r = 90.91 MHz and 
07^/277 = 897.19 MHz, respectively. 

The isotropic chemical shift Siso = — >-^n)/^n is 

a parameter characterizing magnetic shielding effects of 
an atom due to cyclotronic motion of local charges in a 
fixed magnetic field, with w being the actual, measured 
resonance frequency of the nuclei in a given envi¬ 
ronment. This parameter provides detailed information 
on the local chemical environment at the atomic level. 
The temperature dependence of the Siso is shown in 
Fig. 1(a). The signal from PE is shown by black 
squares and that of AFE by red dots. One observes that 
Siso exhibits a continuous behavior in both high and low 
temperature phases but with different slopes, separated 
by a sudden change around the transition temperature 
T]\j. In the narrow transition region one observes the co¬ 
existence of two maxima in the cross polarization and 
magic angle spinning (CPMAS) spectra as detailed in 
Fig. 1(b). Within about 2-3 K of T)v the crystal lattice 
consists of both phases. The PE and AFE phases were 


clearly resolved owing to the exceptionally high spectral 
resolution^® by MAS at 21.1 T (peak widths of about 0.1 
ppm). Note that at any given temperature, the NMR 
peak areas represent the relative abundance of the PE 
and AFE fractions in the crystal lattice. The insets show 
such spectra at Tm = 215.2 K, where the areas of the two 
peaks are about equal, while below and above the transi¬ 
tion temperature the areas are different. It is evident the 
transfer of intensity between the AFE peak (light red) 
and the PE peak (dark gray) as the temperature is low¬ 
ered through Tpf. This shows that one phase is consumed 
by the other, reversibly, with very little temperature hys¬ 
teresis if the temperature variation is slow, about 0.1 K 
per 10 min, as T ^ T^y. 

The high resolution of the peaks allows us to measure 
the spin-lattice relaxation time Ti of the two phases out¬ 
side but also inside the coexisting temperature regime. 
A standard inversion recovery method was used, and the 
signal recovery to the equilibrium following an inversion 
pulse was recorded as a function of recovery time r. Two 
examples are shown in Fig. 2(a), one for each phase at 
Tpf. Briefly after CP, the enhanced ^®N polarization is 
flipped to the —Z axis (that is, antiparallel to the applied 
field). The ^®N magnetization then relaxes back to 

its equilibrium (i.e., +Z axis) along the longitudinal di¬ 
rection. This recovery process can be characterized by 
Ti as follows. 


= 1 - aexp(-T/Ti) (1) 


where a is the cross polarization efficiency. Theoretically 
a could be as large as 10 (~ the ratio of ^H and 

i^N gyromagnetic ratios, respectively), but the exper¬ 
imental setups and relaxation properties could greatly 
affect the CP efficiency^®. In our experiments, a was 
about 4-4.5. As shown in Fig. 2(a), a single exponen¬ 
tial fitting yields different Ti values. Interestingly, for 
the peak representing the PE phase, Ti was found to be 
1.03±0.05 s, while for the AFE peak it was significantly 
shorter, 0.82±0.04 s, both measured at the same temper¬ 
ature (T)v). 

The nuclear relaxation time Ti is an important tool in 
studying phases and their transition. The temperature 
dependence of the ^^N Ti from 295 K down to 180 K is 
shown in Fig. 2(b). A clear discontinuity occurs in the 
Ti data near Tjv ^215 K, implying that the two phases 
have distinctly different motional correlation times. The 
Ti data of Fig. 2b were analyzed using the Bloembergen- 
Purcell-Pound (BPP) model®®: 


1 _ Ctc 3Ctc 

Ti 1 + {UJH - 1 + 

_OCtc_ 

1 -I- {loh + 


( 2 ) 


Here, Tc is the correlation time describing the dynamic 
process and C = ss 4.86 x 10® s"® is 

the dipole-dipole relaxation constant, which is associated 
with the nitrogen-proton distance (N-H bond distance) 
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FIG. 2. (Color online) (a) Normalized magnetization 
of the PE (squares) and AFE (dots) phases vs the in¬ 
version recovery time r at Tn- Single exponential fittings 
(lines) yield Ti of 1.03 and 0.82 s for the PE (black) and AFE 
(red) phases, respectively. The inset is a zoom of data in (b) 
as indicated by the blue rectangle, (b) Temperature depen¬ 
dence of T\ in ADA shown in logarithmic scale for the 
vertical axis and reciprocal scale (l/T) for the temperatnre 
axis, to put in evidence the linear behavior [see Eqs. (2) and 

(3)]. 


Temperature (K) 

FIG. 3. (Color online) Correlation time as a function of tem¬ 
perature, for the AFE (dots) and PE (squares) phases shown 
in logarithmic scale for the vertical axis and reciprocal scale 
(l/T) for the temperature axis. The lines show Arrhenius 
fits done using Eq. 3. The inset shows the same data points 
(no fits) in the same representation, inside the order-disorder 
transition (error bars are smaller than dot sizes). 


rNH = 1.031A taken from Ref. 21. Since Ti is of the 
order of a second, this results in Tc ^ Ins ^ 
consequently Eq. (2) can be reduced to ^ IOCtc. 
For thermally activated diffusion, Tc is assumed to follow 
the Arrhenius behavior, 

Tc = To exp ^ (3) 

where tq is the inverse of a vibrational mode related 
to the phase transition, Ea the activation energy, k the 
Boltzmann constant, and T the temperature. The exper¬ 
imental data in Fig. 2(b) were fitted to the above 
equations, with Ea and tq being treated as variables. The 
slopes of the linear fittings InTi vs T~^, as indicated by 
the solid lines in Fig. 2(b), give different Ea for the PE 
and AFE phases, of 1260 K (or 10.2 kJ/mol) and 780 K 
(or 6.3 kJ/mol), respectively. Similarly, the Arrhenius 
fits given in Fig. 3 give tq equal to 0.3 and 11.6 ps for the 
PE and AFE phases, respectively (and the same values 
for the Ea parameter). The difference in Ea correlates 
well with the fact that the transition induces a signifi¬ 
cant change of the crystalline structure. We ascribe the 
increase in tq as due to the effective “soft” mode lead¬ 
ing to the AFE transition becoming frozen below T/y. 
This indicates that the slowing down of the NH^ tor¬ 
sional motion triggers the phase transition, likely due to 
the formation of new N-H- • • O bonding^^’^®. Moreover, 
as shown below, we can extract microscopic information 
on the size and mass exchange of the ordered phase nan¬ 
oclusters. 

As mentioned earlier on, a major goal of this study 
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FIG. 4. (Golor online) (a) 2D exchange spectrum of ADA 
at 215.2 K with a mixing time of 1 s. The presence of the 
cross peaks indicates a dynamic exchange of the nuclei from 
the AFE clusters to the PE clusters, (b) Schematic of the 
AFE clusters in the PE clustered lattice in going from the PE 
phase (black) to the AFE phase (red). 


was to measure the spatial dimensions of the pretransi- 
tional clusters just above Tjv, their lifetimes and whether 
there is any exchange between the atomic units on the 
periphery of the clustered domain and the PE domains. 
In NMR, such chemical exchange experiments are car¬ 
ried out by introducing a mixing time in between two 
7r/2 pulses, preceded by an evolution time fi, and fol¬ 
lowed by a detection time t 2 - This pulse sequence^^ al¬ 
lows us to record the measured signals during fi _2 (the 
chemical shifts) when the two phases evolve sepa¬ 
rately. During the mixing time, however, the chemical 
shifts start to crosstalk if their corresponding spins (or 
atoms) are involved in a dynamic process, which results 
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in off-diagonal cross peaks after a two-dimensional (2D) 
Fourier transform. Figure 4(a) shows a typical 2D 
chemical exchange spectrum of ADA at 215 K (^ T^r) 
recorded with a mixing time of 1.0 s. The two axes repre¬ 
sent the Fourier transforms of the signals measured dur¬ 
ing the times ti (vertical) and ±2 (horizontal). At this 
temperature, within 0.1 K of T/v, the coexistence of both 
the PE and AFE regions is evidenced by the diagonal 
peaks (at -16.1 and -15.9 ppm), so marked. The 2D spec¬ 
trum also clearly shows the presence of cross peaks, albeit 
small, ~5% as compared to the diagonal peaks. The ex¬ 
istence of the cross peaks implies that the PE and AFE 
clusters are in a dynamic equilibrium such that there is 
a diffusion^^ of molecules between the two phases, and 
moreover that a temperature gradient across the small 
sample is not the cause of the cluster formation. This 
latter was also checked by carrying out measurements on 
a crystal half the size of the previous one, and no signifi¬ 
cant change in the extent of the peak coexistence temper¬ 
ature was observed. We interpret the appearance of the 
cross peaks as showing that some of the PE clusters are 
converted into the AFE domains and vice versa although 
such intercluster exchanges are on the slow time scale of 
the order of the mixing time, ~1 s, at Ttv. Obviously the 
cluster sizes and lifetimes should decrease as the crystal 
temperature is raised well above Tjy. Figure 4(b) shows 
the schematic of the AFE clusters in the PE clustered 
lattice in going from the PE phase to the AFE phase as¬ 
sessed from the relative intensities of the diagonal and 
cross peaks in Fig. 4(a). The cross peaks yield the rel¬ 
ative abundance of the interfacial volume between the 
PE and AFE clusters, i.e., the molecules that are ex¬ 
changing between the PE and AFE domains. These are 
represented by the white boxes in a 2D representation 
shown in Fig. 4. The diagonal peaks yield the relative 
abundance of the PE (black) and AFE (red) domains. 

We then carried out a simple calculation of an estimate 
of the length scale of the clusters as they start to appear 
on the NMR timescale and resolution. Since the AFE 
crystal structure is an orthorhombic (P2i2i2i) phase, 
it is reasonable to assume that the AFE clusters are of 
cubic shape with a side length I, giving a total surface 
area of QP and a cluster volume of If the thickness 
of the cluster surface is assumed to be about the size of 


the unit cell, i.e., ^ 770 pm, the ratio of the boundary 
between the PE and AFE clusters and the AFE cluster 
volume becomes 770pmx6/^/Z^ = 4620 pm/L At Tjv, the 
PE and AFE clusters give about the same signal intensity 
(diagonal peaks), so that the chances of the signals from 
the boundary area going to the PE and AFE clusters are 
about the same. Thus, the ratio of the observed cross 
peak intensity versus the diagonal peak is proportional 
to 2310 pm/L From Fig. 4(a), the cross peak intensities 
are about 5% of the diagonal peak intensities, yielding an 
average cluster size of ~45 ± 5 nm. Therefore, as the PE 
lattice approaches its transition temperature, it develops 
finite size clusters that represent the symmetry of the 
low-temperature phase (i.e., AFE phase), and these AFE 
clusters increase in size and numbers to fill the entire 
lattice as the temperature goes down below 215.2 K, as 
illustrated schematically in Fig. 4(b). 

Concluding, we detail the structural and dynamic as¬ 
pects of this antiferroelectric transition and give evidence 
for the co-existence of two phases during such an order- 
disorder transition. This was achieved using the excep¬ 
tionally high spectral resolution provided by MAS NMR 
at 21.1 T. The quantitative estimate of the AFE clusters 
size is ~45 nm and their lifetimes are ^1 s at T/v-1-0.1 K. 
These quantitative data should serve for the develop¬ 
ment of precise models of AFE phase transitions. The 
reported NMR technique should also open a window on 
other phase transitions, especially for all nonconducting 
and nonmagnetic solids for which high-resolution solid 
NMR is a sensitive structural and molecular dynamics 
probe. 
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